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Abstract 

Long-range correlations between multipUcities in different rapidity windows in hadron- 
nucleus collisions are analyzed. After recalling the standard results in the probabilistic 
model, we study them in the framework of perturbative QCD. Considering interacting BFKL 
pomerons in the form of fan diagrams coupled to a dilute projectile, analytic estimates are 
done for very large rapidities. The correlation strength results weakly depending on energy 
and centrality or nuclear size, and generically greater than unity. Finally, we turn to the 
Color Glass Condensate framework. For a saturated projectile and considering the most 
feasible experimental situation of forward and backward rapidity windows symmetric around 
the center-of-mass, the resulting correlation strength turns out to be larger than unity and 
shows a non-monotonic behavior with increasing energy, first increasing and then decreasing 
to a limiting value. Its behavior with increasing centrality or nuclear size depends on the 
considered rapidity windows. 



1 Introduction 

Long-range correlations have been attracting much attention since long ago in the region 
of both low and high transverse momenta of secondaries. At low momenta the color string 
picture [1] with fusion and percolation effects [2] has been extensively applied [31 H]. In the 
semihard region the Color Glass Condensate (CGC) picture has lately been used OE]. One 
can separate long-range correlations into a contribution from purely hadronic collisions and 
another coming from the effects generated by the heavy nucleus target or/and projectile. 
Obviously the first contribution can hardly be treated in a more or less rigorous theoretical 
framework due to the essentially non-perturbative structure of the hadron. Heavy nucleus 
collisions, on the other hand, present more opportunities in this sense, due to their compar- 
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atively simple structure in terms of constituent nucleons. Single inclusive cross-sections with 
participation of nuclei can be easily found even in the framework of interacting pomerons, 
both of the old local type and of the sophisticated Balitsky-Fadin-Kuraev-Lipatov (BFKL) 
type. In the limit of very heavy nuclei they are expressed in terms of the sum of the corre- 
sponding pomeron fan diagrams. However long-range correlations require also knowledge of 
the double inclusive cross-section, for which the situation is more complicated. In the purely 
eikonal approach they can also be easily calculated from the known single and double ele- 
mentary (hadron-nucleon) inclusive cross-sections. However with interacting pomerons this is 
only possible for the hadron-nucleus case. Double inclusive cross-section for nucleus-nucleus 
scattering mediated by interacting pomerons involves a complicated set of diagrams, exact 
summation of which does not look feasible. So hadron-nucleus collisions present a subject 
better suited for theoretical discussion of long-range correlations in the nuclear background. 

In this paper we study long-range correlations in hadron-nucleus collisions in the hard 
domain. We shall use two different approaches to this problem, which treat different range 
of energies. At very high (asymptotic) energies we shall rely on the perturbative Quan- 
tum Chromodynamics (QCD), which predicts that the interaction is mediated by interacting 
BFKL pomerons while treating the hadron as a dilute object through the whole evolution. 
At smaller energies, when evolution of the gluon density is not complete, we shall apply the 
Color Glass Condensate approach, in which the fast nucleus is represented by a strong classi- 
cal gluon field [7] . This approach has been lately used for a qualitative study of correlations 
in nucleus- nucleus collisions with promising results [6j. 

The paper is organized as follows. In the next section, to have a basis for comparison 
and discussion, we shall recall some basic predictions for the long-range correlations in hA 
collisions, which follow from the straightforward probabilistic approach. Section 3 is ded- 
icated to the perturbative QCD approach at very high energies. Section 4 presents some 
qualitative predictions from the Color Glass Condensate. Finally we discuss our results in 
the Conclusions. 

2 Probabilistic treatment 

In this section we shall study long-range correlations in hA collisions, as they follow from 
purely probabilistic considerations. This approach is realized in the Glauber-eikonal descrip- 
tion of hA interactions and also reproduced in Regge approach with non-interacting pomerons. 
Our basic instrument will be the standard Glauber expression for the cross-section cr„ for n 
inelastic interactions of the projectile inside the nucleus: 




(1) 
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Here a is the elementary (hadron-nucleon) cross-section, and T{h) is the standard nuclear 
profile function normalized to unity. From this expression wc immediately derive expressions 
for the single and double inclusive cross-sections for hA collisions, Ji and J2 respectively. At 
fixed impact parameter h 

Jiiy,k)=An{y,k)T{b) (2) 

and 

J2{yi,ki;y2,k2) = Aj2iyi,ki;y2,k2)T{b) + A{A - l)ji{yi,ki)ji{y2,k2)T\b). (3) 

Here and in the following we denote by small letters the quantities which refer to the el- 
ementary hN collision. So j'l and j2 are the single and double inclusive cross-sections for 
hadron-nucleon collisions; y and k denote the rapidity and transverse momentum of the 
produced particle. 

To pass to the corresponding multiplicities we have to integrate over b and divide by 
the total inelastic hA cross-section S*". To do this we have to choose a form of the profile 
function T{b). We use the simplest choice of the constant nuclear density within a sphere of 
radius Ra = A^/^Rq, which gives 



T{b) = ^^-A , (4) 



where Va is the nuclear volume. With this profile function we find 
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For large A, with a good precision, the inelastic hA cross-section S'" = TrR\, so that we 
finally find the single and double hA multiplicities. Mi and M2, respectively as 

Mi{y,k) = A^/^mi{y,k) (6) 

and 

9 

M2{yi,ki;y2,k2) = A^'^m2{yi,ki;y2,k2) + -A'^'^ni{yi,ki)mi{y2,k2). (7) 

o 

With these expressions we can pass to correlations. The strength of long-range correlations 
in nuclear collisions is standardly determined by the coefficient 

_ {NfNb) - {Nf){Nb) 

where Np and Nb are the numbers of particles produced in two rapidity windows, separated 
by a sufficiently large rapidity interval ('forward' and 'backward'). Note that in the asymmet- 
ric hadron-nucleus case, there is another correlation strength defined with {Np) — {Np}^ in 
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the denominator. As to the particle transverse momenta, they may be both taken integrated 
over the whole phase space or restricted to specific parts of it (even practically fixed). This 
circumstance plays no role for the following derivation. The average numbers which figure in 
([8]) are expressed via the multiplicities as follows 

{NFiB)) = I dT^^^^M,iy,k), (9) 

{NfNf^b)) = J dT[dT^^''^M2{yi,h;y2,k2), (10) 

where dr^ and dr* denote integration over y and k in the forward and backward windows. 

Similar quantities for the elementary hN collisions will be denoted by small letters. So 
for hN collisions the correlation coefficient is determined as 

_ {npUB) - {nF){nB) , . 

{nl)-{np)^ ' ^''^ 

where the averages are defined as in ([9]) and ([TO]) with multiplicities and /i2- Using 
relations ^ and d?]) we can express averages for the nucleus target via the same quantities 
on the nucleon target to obtain: 

{nFnB) + lA^''^{nF){nB) 

If we define the dispersion squared in the forward window for the elementary collisions as 

d^ = {nl) - {npf, (13) 

then we find 

B = ^ ^ (14) 

For symmetric windows (relative to hN collisions) we have (np) = {ub) so that ()14p simplifies 
to 

Q = ' (15) 

d2+ (1^1/3 + i)(n^)2 ■ 

If the dispersion squared is much smaller than the particle number squared, which is expected 
for large enough energies and windows, we find an approximate expression 

As a result we find that in the theoretical limit ^ ^ oo the nuclear coefficient B tends to 
unity, the value of the elementary coefficient h having no importance. So, in this limit, long- 
range correlations are exclusively a consequence of the nuclear environment. Of course for 
realistic nuclei the term (1/8)^4^/3 ^ot large but smaller than unity so that the coefficient B 
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results noticeably smaller than unity. However we also see that its value only weakly depends 
on A and is mainly determined by the relative dispersion d/ {np) in elementary collisions. All 
these effects are due to the presence of the second term in the nuclear multiplicity M2, Eq. 
([7j). Note that the limiting case yl = 1 is achieved through the substitution A"^/^/8 ^ — 1 
in (114p . so that the second terms in both the numerator and denominator of this equation 
vanish and B passes into b. 

3 Long-range correlations in high-energy hA colhsions in the 
perturbative QCD 

3.1 Formalism 

In the perturbative QCD at high energies the interaction between the incoming hadron and 
the nuclear target is realized by an exchange of BFKL pomerons, which interact between 
themselves via the triple pomeron vertex. For a heavy nucleus target with A ^ 1 and a 
point-like projectile this interaction is described by a set of pomeron fan diagrams, which are 
summed by the non- linear Balitsky-Kovchegov equation [HI O [10] . If we denote this sum for 
a fixed impact parameter 6, zero total gluon momentum and intergluon transverse distance 
r as r, b) then function (f){y, r, b) = ^{y, r, b) /(27rr^), transformed into momentum space, 
satisfies the equation 

riA^in n hi] 

= -H4>{y,q,b)-^^{y,q,b), (17) 



dy 

where H is the BFKL Hamiltonian (see e.g [11]) and y = ay is a rescaled rapidity with the 
standard notation a = OsNc/tt. The physical meaning of (p is provided by its relation to the 
gluon density in the nucleus: 

d(xG{x,k^,b)) n2 I 

where 

h{y,k,b)=k'vl<P{y,k,b) (19) 

and y = ln(a;o/x) with xq some initial value of x, usually taken ~ 0.01. Function h satisfies 
a normalization condition 

'^^^ h{y,k,b) = l. (20) 



27rA:2 

Numerical calculations [12] show that starting from y ~ 2 it develops a scaling structure 



%,A:,6) = /i(e) = 0.295 exp(-CV3.476), i = \n[k/Qs{y,b)) , (21) 
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where Qs{y,b) is the so-called saturation momentum, which grows exponentially with y: 



Ai = (2.0 H- 2.3)a and a is a numerical constantly. 

An immediate physical application of this framework is to deep inelastic scattering (DIS), 
with a highly virtual photon as a projectile. A hadronic analogue of this may be the 'onium', 
that is a quark-antiquark system of a very short dimension. Realistic hadronic projectiles 
are not point-like and do not allow for the perturbative treatment. So application of the fan 
diagram approach to their interaction is strictly speaking not very well justified. For this 
reason in the following we have to keep in mind an approximate character of our derivation, 
which assumes that, as in DIS, the projectile hadron interacts with the pomeron only once. 
Under this approximation the total hadron-nucleus cross-section is given by 



where p{r) is the color density of the projectile and Y is the overall rapidity. It may be 
illustrated graphically as shown in Fig. 1, where the circle with the attached line correspond 
to function <I>. For a normalizable density p{r) and finite nucleus (with the profile function 
dH)) cross-section S = S*" turns out to be purely geometric: S = 2t:R\. 



3.2 Single inclusive cross-sections 

The single inclusive cross-section is given by the two diagrams shown in Fig. 2 a and b. They 
correspond to production of the observed particle either from the pomeron which couples to 
^This dependence on impact parameter and nuclear size was obtained [T2] for a realistic profile function. 
Also the value of Ai depends on the rapidity window of the fit, and is slightly smaller than the asymptotic 
theoretical expectation, Ai = 2.44a. See detailed discussions of these aspects in e.g. [13) where a dependence 
oc A^^^ is obtained for a cylindrical nucleus. None of these considerations alter the conclusions obtained in 






(23) 




Figure 1: Diagram contributing to the total hA cross-section. 



this Section. 
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the projectile or from the vertex of its sphtting into a pair of pomerons \14:\ [T5j . Ah other 
possibihties are canceled by the Abramovsky-Gribov-Kancheli cutting rules [16]. The sum of 





Figure 2: Diagrams contributing to the single inclusive hA cross-section. 



these two contributions leads to the single inclusive cross-section at fixed b [17 

Ana 



Ji{y,k,b) 



fc2 



■ J (fre'^''V'^P{Y - y, r)V^ (2$(y, r, b) - ^'^{y, r, 6)) , 



where P{Y — y, r) is the pomeron coupled to the projectile: 

P{Y -y,r) = J d\'p{r')G{Y - y, r', 



(24) 



(25) 



with the nucleus at y = and the dilute projectile at y = Y , and G is the BFKL forward 
Green function [llj. Performing the differentiations and passing to momentum space one 
obtains 



Ji{y,k,b) 



Ana 



cPqh^'^^iY — y, q)w{y, k — q,b). 



(26) 



Here h^^^ (y, q) is a function analogous to h{y, q, r) for the pomeron, that is a Fourier transform 
of V2p(y,r)/(27r). Function w{y,q,b) is defined via h{y,q,b): 



^ /■,2 h{y^Qi,b)Hy^Q - Qi^b) 



(27) 



w{y,q,b) = ^ I d^qi- ^ 

It satisfies the same normalization condition (j20p and has the same scaling property (j2ip . 
although with a shifted maximum and slope in its ^ dependence: 

w{y, k, b) = 0.358 exp ( - 0.402(^ - 0.756)^) . (28) 

Of course the shift in the maximum corresponds to a universal enhancement of the value of 
the saturation momentum. 

Numerical calculations of the single inclusive cross-section along these formulas were per- 
formed in [T7|. However such calculations for the double inclusive cross-sections look very 
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difficult, so that we shall instead use analytical estimates obtained in the asymptotic regime 
when both Y — y and y are large. Then we can use the well known asymptotic expressions 
for the BFKL Green function to obtain an explicit expression for function h^^\ Repeating 
the calculations done in [14J we obtain 



8a 



TT 



F{y,k,b), 



(29) 



where A = 41n 2a is the BFKL intercept, /? = 14C(3)a, Rp is the radius of the projectile and 

(30) 



- — w{y,k - q,b) = Qs{y,b)f{k), k = — — -. 
2TTq Qs{y,b) 

Function f{k) obviously reduces to a constant /(O) when k <^ Qs- In the opposite case when 
k ^ Qs it behaves as X/k. Numerical calculations give values for f{k) illustrated in Fig. 3 
with 

/(0) = 3.97, A = 54.6. (31) 
So for momenta considerably below the saturation momentum we find 




Figure 3: Function /(e^). 



My,k,b) = ^Rpf(0)e^iy-y)^-^^^^Q,(y,b) (32) 
and for momenta considerably above the saturation momentum 



Ji{y,k,b) = ^RpXe^(^-y^^^^^^^^Q\y,b). (33) 

The found cross-section grows exponentially with the overall rapidity Y, which just reflects 
the growth of the pomeron directly coupled to the projectile. One expects that for an extended 
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projectile this growth would be finally tempered when more than one pomeron are coupled 
to it (see Fig. 4 ), contributions which are damped by powers of the small coupling constant 
for a point-like projectile. As we shall see, in the expression for the correlation coefficient the 
growing factors cancel, so that the resulting coefficient depends on Y only weakly. For this 
reason we may hope that our formulas remain applicable also for realistic hadrons. Another 




Figure 4: Diagrams not taken into account for the total cross-section for a point-like projectile. 

interesting feature is the peculiar dependence on the nuclear factor (^AT{b)^ , with p = 2/3 
for k <^ Qs and p = 4/3 for k ^ Qg, which is of primary importance for the correlations. 
Integration over b and division by S leads to multiplicities in these two regions of k. For 
k<.Qs 

M,iy,,)^A'''y^R,me^y-.>^ZJIZ^ (34) 

where 

/ 9 \i/3 

7 = «(i^) (35) 



and e = A — Ai, and for k ^ Qs 



with ei = A - 2Ai. 

3.3 Double inclusive cross-section and correlations 

The double inclusive cross-section is described by a set of diagrams shown in Fig. 5 (a few 
evident additional diagrams are not shown). They are quite complicated, especially since 
the cut vertex appearing in the diagram Fig. 5 / is different from the uncut one |18j . As 
mentioned in the Introduction, a detailed calculation of all the contributions does not seem 
very realistic. However at high values of all rapidity distances, Y — yi,Y — y2,yi,y2 ^ 1, of 
all the contributions the dominant ones correspond to Figs. 5 c-e, in which the upper vertex 
can have rapidities up to Y, so that the two pomeron lines below give the maximally growing 
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exponential factor exp (^A(2Y — yi — 2/2)) • The relative weights of all other contributions is 
exponentially damped. The study of contributions from all the diagrams of Fig. 5 c-e can be 




a 






Figure 5: Diagrams contributing to the double inclusive cross-section in perturbative QCD. 

done by the same method which was used in [14] for the diagram of Fig. 5 c. Contribution of 
the diagrams 5. d and e can be taken into account by changing function h{y, q, r) to w{y, q, h) 
in all the formulas. In this way we obtain the double inclusive cross-section as 



Myi,h;y2,h,b) 



4 



:{r )pe 



IT 



-Fiyi,ki,b)Fiy2,k2,b), 



(37) 



where F is the same function (|30p. Further simplifications, similarly to the case of sin- 
gle inclusive cross-sections, can be made when ki and k2 are either much smaller than the 
saturation momentum or much larger than it. In the case when both ki ^ Qg and k2 <C Qs 

4 a2 ' 



'^^ In 2 kjk^ 



(r2)p/2(0)e^(2y-y,-j;,) 



TT 



/3(2y - yi - y2) 



Qsiyi,b)Qsiy2,b). (38) 



If both ki ^ Qs and k2 ^ Qs then 



J2 



4 a' 



IT 



In 2 kfk^ 

Finally in the case ki » Qs and k2 ^ Qs we find 
4 a2 



Rf'yv ^Q'siyi,b)Qliy2,b). (39) 

P{2Y - 2/1 - 2/2) 



J2 



TT 



-Qliyi,b)Qs{y2,b). (40) 



ln2kfkf '^"^--'^ Y /3(2y _ _ ys) 

Integration over b and division by S leads to the corresponding multiplicities Al2{yi, ^i; 2/2) ^2)- 
In the case when both ki <^ Qs and k2 <^ Qs 



4 



10 ' ln2 kfk^ 



{r^)pf\^)e^^''-<y^+y^\ 



vr 



Pyiy2{2Y - yi - ^2) 



(41) 
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If both ki » Qs and k2 ^ Qs then 



Finally in the case ki » and k2 ^ Qs we find 

^2 = A'/'j^^{r')pXfiOW'e''^'^-^^y^-^y^ I , / (43) 
4 In 2 A;Jfc| y [3y{y2{2Y - yi - y2) 

It is convenient to introduce a ratio 

, , . M2{yi,ki;y2,k2) 
Mi(yi,fei)Mi(?/2,fc2) 

for which for all three limiting cases considered above we obtain a simple expression: 



mi. y\ r ^ ^"^'^P P{Y-y,){Y-y2) 

R{yi,h]y2,k2) = C-———^J — — — , 45 

16 In 2 Rj, y TT{2Y-yi-y2) 

where for the cases ki,k2 ^ Qs, ki,k2 ^ Qs and ki » Qs,k2 ^ Qs the coefficient C is 
32/15, 50/21 and 20/9 respectively. If the two produced jets of hadrons are taken symmetric 
in the cm. system for hN collisions with the rapidity distance y then 

yi = \{Y + y), y2 = liY-y), (46) 

and 

R{yi, ki;y2, k2) = C-^^J ^ r(yi, ya), (47) 
62 in 2 Hp V 

with 

/ y2 



r(yi,2/2) = y^-|7, r(yi,yi) = yy^, T{y2,y2) = VYT^. (48) 

Eqs. (|15]) . (d?!) and (|18]) are the central result of this Section. 
In terms of the ratio R the correlation coefficient is given by 

^ ^ MijyB, ks) Rjyp, kp; ys, kn) - 1 .^g. 
Mi{yp,kp) R{yp,kp;yp,kp) - 1' 

The coefficient in Eq. (j47p is very small e.g. ~ 0.1 for ~ 0.2. Thus both the numerator 
and denominator in this equation are negative except for very large energies (Y > 100 in 
the mentioned example). Furthermore, R may be smaller than one for some small window 
in y and Y, but it is generically larger than one. To illustrate it, at large Y and fixed y we 
conclude from (1470 that R is independent of y in all cases. So if kp and kp have the same 
order of magnitude (either much smaller or much larger than Qs) the second ratio in ()49|) is 
equal to unity. In this case we have a simple result 

B = ^^P44- (50) 
Miiyp,kp) ^ ' 
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For the most important case from the practical point of view, kp, ks <C Qs, we then find 

B = ^e-, (51) 

or for windows symmetric also in the phase volume of transverse momenta simply 

B = e'y . (52) 

The concrete value of e depends on the chosen value for Ug- With a typical value Ug = 0.2 
we find e ~ 0.1 0.15. 

It is not difficult to obtain predictions for B for all other theoretically possible cases. If 
both kp and ks are much larger than the saturation momentum we obtain the correlation 
coefficient (|52p with e substituted by ei. In still more exotic cases, when one of the momenta 
is much smaller and the other much larger than Qg, the second ratio in (j49p begins to depend 
on Y non-trivially because of different coefficients C in (j47p . Also a nontrivial dependence 
on A appears, due to different powers of A in (j34p and (|36p . The explicit formulas can be 
easily written using (I34p . (j36p and (j47p . We do not present them due to a small probability 
of the corresponding experimental setup. 

As we see in all cases the correlation coefficient turns out to be different from the proba- 
bilistic predictions. For symmetric windows it is independent of j4, generically greater than 
unity and grows (rather slowly) with the rapidity distance. 



4 Color Glass Condensate 

In this Section we follow the lines in [5l [6] to obtain an expression for the correlation strength 
B in hadron- nucleus, pA collisions, considering the hadron as a saturated object with some 
saturation scale Ql^p{y) < Ql a^v) ^ ^ > 0, as done in [191 EQ]. For the multiplicity 
of produced gluons, one gets in a small overlap area (with a <^ Rq corresponding to the 
correlation length of the classical fields) between projectile and target: 



/dN\ Qs^mniy) Qlrmniy) = { <3s,p(y) ' <5s,A(y) } • (53) 

\ dy / as {Qs,min{y)) L '^^ ' J 

After integration over impact parameter one gets an overlap area S i.e. 

\ dy / as {Qs,rriin{y)) 

such that SQ'l ^ oc Npart,hi the number of nucleons from hadron h participating in the collision 
{Npart,p = 1), see also [20]. 

The numerator in Eq. Q (coming from diagrams like that in Fig. 6 a) results in 

dN dN\ /dN\/dN\^ Qlminiyp) QlrmniyB) ^4 ^^^^ 



dypdyB/ \dyF/\dyB/ as {Qs,min{yF)) as {Qs,min{yB)) 



a 
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Figure 6: Diagrams contributing to the double inclusive cross-section in the CGC. The black 
dot and the cross correspond to the projectile and the target from which the classical fields, 
represented by the solid lines, come. The dashed lines represent the emitted gluons. 

Now, for the integration over impact parameter we consider that a} ~ \IQ\^^^{yF-,yB\ 
QlmaxiyF,yB) = max | , ^ (?/i7') , ^(ys) , ^ (t/s) | , which results in 



dN dN 



dN\ / dN 



S 



(56) 



.dypdyB/ Xdyp/Xdys/ Os {Qs,min{yF)) Os {Qs,min{yB)) Q%maxiyF,yB) 

For the denominator in Eq. ([8]) there is an additional piece coming from diagrams like that 
in Fig. 6 b, 

\/dNdN\ /dN\/dN\]' ^ , , , , „^ 

As discussed in [6J, this piece is ©(a^) suppressed compared to ([56l) . It contains an expo- 
nential damping factor, with k ~ 1, which motivates its absence in the numerator for large 
enough yp — ys- Also note that for symmetric AA collisions, both (j56p and ()57p are ©(A"^/^) 
suppressed compared to the square of ([54ll . as discussed in [5]. 

Neglecting any possible interference between these two kinds of diagrams and for large 
enough yp — ys, the final expression for Eq. ([8]) readj^ 



B 



as iQs,min{yF)) Ql,max{y F ,yF)Qlmin{yB ) 



l + Cal {Qs,min{yF)) 



Qimax{yF,yF) 
Qlminiyp) 



as {Qs,min{yB)) QlmaxiyF , yB)Ql,niniyF) 

' (58) 

where c is a constant. Note that for AA collisions and symmetric intervals yB = —yp, this 
expression reduces to that found in [6], = [1 + ca^] ^ , from which the correlation strength 



An error in this formula in a previous version of this manuscript, led to wrong conclusions in this Section. 
We warmly thank J.G.Milhano for pointing it out. 
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was argued to increase with centrahty and energy of the colhsion. For the asymmetric case 
Eq. (f58]) . these behaviors depend on the considered rapidities yF,yB- 
To iUustrate this, let us take the parametrizations [20] 

with Qo a constant with dimension of momentum squared, Ecm the collision energy in the 
center-of-mass system, A ~ 0.3, and rapidities defined in the center-of-mass system in which 
the projectile is located at y = Y/2 and the nuclear target at y = —Y/2. The solution of the 
equation 

Qlpivc) = QUivc) (60) 

defines [20j a critical rapidity yc — —4 — 3 < such that for y < yc, Qlminiv) — Ql Aiv)i 
while for y > y^ Q'iminiy) — Ql,p{y)- Let us examine several situations: 

• For yc < Vb = —yp - the most feasible situation from the experimental point of view, 

Qs,min{yF) = Qs,p{yF), Qs,min{yB) = Qs,p{yB) and Q s ,max{y F , y b) = Qs^iyp), the 

correlation strength B is generically greater than 1 and decreases with increasing A. 
It shows a non-monotonic behavior with increasing energy, first increasing and then 
decreasing towards a limiting value e'^'^^^ > 1 independent of A. This case coincides 
with the one in ()52|) (with e A and y ^yp), which can be easily understood as 
the consideration of a dilute projectile in the previous Section is equivalent to the limit 
yc — > — oo here. 

• For yB = -yp < yc, Qs,min{yF) = Qs,p{yF), Qs,min{yB) = Qs,A{yB) and Qs,max{yF,yB) 

= Qs^A^yF), the correlation strength B is generically greater than 1 and increases with 
increasing A. Again, it shows a non-monotonic behavior with increasing energy, first 
increasing and then decreasing towards a limiting value A^'^~^'^/^ > 1 independent of 
yF = -yB- 

• For yB < yF < yc situation within the kinematical reach of the Large Hadron Col- 
lider (LHC) with e.g. ys = -7 and yp = -5), Qs,min{yF) = Qs,AiyF), Qs,min{yB) = 
Qs,AiyB) and Qs,maxiyF,yB) = QsAys), the correlation strength B is generically 
smaller than 1 and increases with increasing A. Opposite to the two previous cases, it 
shows a monotonic increase with increasing energy towards a limiting value e"'^'**'-^^"^^^ 
< 1 independent of A. 

Let us note that these estimates have been done for truly asymptotic values of energy and 
nuclear size and for a concrete choice of the parametrization for the saturation scale, the 
behavior for the experimentally accessible situation depending on this concrete choice and 
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on the value of parameter c. Nevertheless, the fact that the correlation strength is larger or 
smaller than one is generic. 

The results in this Section should be applicable for transverse momenta of the order or 
smaller than the corresponding saturation scales. In any case, the correlation strength is 
larger than 1 for the most feasible experimental situation of forward and backward rapidity 
windows symmetric around the center-of-mass, a behavior which coincides to that generically 
found in the previous Section, see e.g. (f52|l . 

5 Conclusions 

In this paper we analyze the long-range rapidity correlations in hadron-nucleus collisions. 
First we recall the standard results in the probabilistic model. The correlation strength B 
shows an increasing behavior with increasing nuclear size or centrality, tending to unity for 
^ ^ oo. 

Next we turn to perturbative QCD. We consider interacting BFKL pomerons in the form 
of fan diagrams coupled to a dilute projectile. After examining the required single and double 
inclusive cross-sections \17\ [T8] , analytic estimates are done for very large rapidities due to 
the difficulties for a complete consideration of the double inclusive density. The correlation 
strength results weakly depending on energy and centrality or nuclear size, and generically 
greater than unity. Note that these results are rigorously applicable to DIS, but require 
certain caution if applied to hA scattering with ordinary hadrons. 

Finally, we turn to the Color Glass Condensate framework. Taking the projectile as a 
saturated object characterized by a saturation scale smaller - as expected - than that of the 
nucleus, we extend to asymmetric collisions the analysis done in [5l [6] for nucleus-nucleus 
collisions. For the most feasible experimental situation of forward and backward rapidity 
windows symmetric around the center-of-mass, the resulting correlation strength turns out 
to be larger than unity and shows a non-monotonic behavior with increasing energy, first 
increasing and then decreasing to a limiting value. Its behavior with increasing centrality or 
nuclear size depends on the considered rapidity windows. 

A note of caution is in order. The correlations considered in this paper are those coming 
from particle production in the initial stage of the collision. Subsequent stages may modify 
the predicted behaviors. In any case, hadron-nucleus collisions should offer a more reliable 
setup than nucleus-nucleus in this respect, as the production of a dense, eventually thermal- 
ized medium is not expected. Besides, hadronic rescattering of the produced secondaries is 
expected to play a little role except in the region close to the rapidity of the nucleus. With 
all these caveats in mind, phenomenological applications of our results to dAu collisions at 
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the Relativistic Heavy Ion Collider and pA collisions at the LHC are the obvious extension 
of this work. 
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